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Abstract. We establish sharp bounds for simultaneous local rotation 
and Holder-distortion of planar quasiconformal maps. In addition, we 
give sharp estimates for the corresponding joint quasiconformal mul- 
tifractal spectrum, based on new estimates for Burkholder functionals 
with complex parameters. As a consequence, we obtain optimal rotation 
estimates also for bi-Lipschitz maps. 



1. Introduction 

A deformation /: M? — > M 2 is called L-bilipschitz if it distorts Euclidean 
distances by at most a fixed factor L 1, 

\x - y\ ^ \f{x) - f(y)\ ^ L\x-y\ for x,y G M 2 . 

Such a map changes length insignificantly, nevertheless it may change local 
geometry by creating (logarithmic) spirals out of line segments. As a simple 
model example consider the logarithmic spiral map 

f(z) = z\zr, zeC = R 2 , 7 eE\{0}. (l.i) 

On the other hand, the constant L imposes constraints on the speed of 



spiralling: for any L-bilipschitz map, see Proposition 6.1 below, the infini- 
tesimal rate of rotation at a point z£K 2 , 

, . r ^g[f(z + t)-f(z)] 

7 (z) =jf(z) = km sup = , 1.2 

t-+o log t 

satisfies \^\ ^ L — \ with equality for (1.1) at z = 0. 
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Note that a map can have non-trivial (i.e. 7 / 0) rotation only at points 
of non-differentiability, hence for bilipschitz maps only in a set of measure 
zero. This leads one to study the rotational multifractal spectrum of these 
mappings, that is, to ask what is the maximal Hausdorff dimension of a set 
where an L-bilipschitz deformation rotates every point with a given rate 

7- 

The following result gives a complete answer, describing the universal 
bounds for the class of bilipschitz mappings. Indeed, these interpolate lin- 
early between the bounds valid pointwise and those valid almost everywhere. 

Theorem 1.1. Suppose f : M 2 — > M 2 is L-bilipschitz and 7 is a real number 
such that \j\ ^ L — Then 

2L 

dim H {z G R 2 : V (z) = 7} ^ 2 - (1-3) 



Moreover, for every such 7 there exists an L-bilipschitz map f : 



p2 



for which the equality holds in (1.3). 



A key observation towards establishing rotational bounds, such as Theo- 
rem [TTJ is that it is necessary to control the higher integrability of complex 
powers of the derivative f z = df jdz. Such an integrability is most naturally 
studied through a holomorphic flow of the map /. 

However, holomorphic flows do not keep the maps bilipschitz. Any orien- 
tation preserving L-bilipschitz / : M 2 — > R 2 solves the Beltrami equation 

L 2 - 1 

f- z = u(z)f z , with |^)| <^— _ a.e. x£R 2 . (1.4) 

Replacing \x by \x\ , a coefficient depending holomorphically on a complex pa- 
rameter A, the equations provide us with homeomorphic solutions which vary 
holomorphically with A - holomorphic motions of Mane, Sad and Sullivan 
|17j . However, in general these homeomorphisms are only quasisymmetric 



(see ( |2.2[ ) below) even if the initial map in (1.4) is bilipschitz. 

Surprisingly, though we are forced to move outside the bilipschitz world, 
this approach leads to a complete description of the multifractal rotation 
spectra. 
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The correct setup for integral estimates of complex powers of f z is within 



the solutions to an arbitrary Beltrami equation (1.4), thus in terms of qua- 
siconformal mappings. Here the result takes the following form (for the 
appropriate concepts see Section [2]). 

Theorem 1.2. Suppose f is a K -quasiconformal map on a domain Q C C. 
Then for any exponent f3 £ C in the critical ellipse 

|/3| + |/3-2|<2-^ (1.5) 



we have 



ft 



e Lj oc (n). 



This result is sharp in a strong sense; it fails for any complex exponent f3 
on the boundary of the critical ellipse as well as outside, see Section |4| 

Theorem |1.2| includes a number of interesting special cases. Taking /3 
real-valued we recover the optimal higher integrability of the gradient of 
a quasiconformal mapping [lj. Other values of (3 lead to new phenomena: 
among them is the optimal exponential integrability of the argument of f z . 

Corollary 1.3. Suppose f is a K -quasiconformal map on a domain Cl C C. 
Then 

AK 

e b|arg/ 2 | g ^ for all positive b < 



K 2 -l' 

Again the integrability fails for some i^-quasiconformal mapping in case 
b = j^zi ■ Moreover, we have similar precise bounds on exponential inte- 
grability for bilipschitz mappings as well, see Theorem |6.2| below. 

To understand the mechanisms of how the integrability breaks down at 



the borderline case in Theorem 1.2, we analyze the situation in a weighted 
setting. This reveals connections with Burkholder functionals and raises 
new questions regarding quasiconvexity. The novelty here lies in creating a 
family of Burkholder- type functionals depending on a complex parameter, 



for which we establish partial quasiconcavity, see Theorem 4.3. In the course 
of doing so we bring into play Lebesgue spaces with complex exponents to 
advance in this setting the interpolation technique of jHj. For details see 
Section HI 

The above approach allows a natural extension to multifractal properties 
of quasiconformal mappings. Even further, in view of |17j . [22j one may 
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consider general holomorphic motions of subsets ficC, that is, maps ^/ : 
Dxi?-}C injective in the z-variable and holomorphic in the A-variable, 
with ^(O, z) = z at the "time" A = 0. Now, however, one needs to take into 



account also the stretching in a manner to be discussed in Section 3.2 We 
arrive at general bounds for the joint rotational and stretching multifractal 
spectra. 

Theorem 1.4. Suppose $:Dxi?->Cisa holomorphic motion of a set 
EcC and that a > and 7 £ 1 are given. 

If X £ D, assume that at every point z G E we have scales rj — > along 
which ^\(z) = ty(A,z) stretches with exponent a, 

]im **\*xiz + r j )-9 x (z)\ = 
j^oo log rj 

and simultaneously rotates with rate 7, 

lh arg(^ + ,)-f A (,)) 
j-+oo log\^ x {z + rj )-^x{z)\ 

Then 



dim(E) ^ 1 + a - —^{l-a) 2 + (1 - |A| 2 )a 2 7 2 - (1-6) 
1^1 

Moreover, for each A £ D, a > and 7GI such that the right hand side 



of ( 1.6) is nonnegative, there exists a set E C C and its holomorphic motion 



for which we have equality in (1.6). 

Holomorphic motions and the study of their geometric properties arise 
naturally in various questions in complex dynamics. It is clear that detailed 
combinatorial or topological information about specific dynamical systems, 



combined with the methods of Theorem 1.4, will improve the bounds above. 

The proofs of the statements in this introduction can be found in the 
text as follows. Theorem \\.2\ and Corollary 1.3 are proved in Section |4j 



Theorem 1.4 and further quasiconformal multifractal spectra estimates are 
treated in Section [5j Finally, Theorem |1.1| is obtained by developing the 
bilipschitz theory in Section [6} In Section [3] we build up the basic framework 
for discussing rotational phenomena through rigorous definitions of various 
branches of the logarithms involved. Section [2] contains prerequisites on 
bilipschitz and quasiconformal mappings. 
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The above and further results in the present work as well as of other 
authors reflect a close interaction between rotational phenomena for bilip- 
schitz mappings and stretching phenomena for quasiconformal mappings. 
This may be summarized as the dictionary of Table [T] between the two, see 
Section [6] for a discussion. 



K-quasiconformal stretching 


L-bilipschitz rotation 


Grotzsch problem 


John's problem 


radial stretching 


logarithmic spiral map 


z\z\ a ~ x 


z \z * 7 


Holder exponent 


rate of spiralling 


1/K ^a^K 


W^L-l/L 


log J(z, /) G BMO 


a vgf z eBMO 


higher integrability 


exponential integrability 


J fc vv loc ' P ^ K-l 


exp(6|arg/*|) G Lj oc , b < £\ 


multifractal spectrum 


multifractal spectrum 


dim^{z : a(z) = a} ^ 1 + a — 


dim H {w : j(w) = 7} ^ 2 (7! 


factoring the radial stretch map 


factoring the logarithmic spiral map 


along a geodesic 


along a horocycle 



Table 1. Quasiconformal stretching versus Bilipschitz rotation 



2. Prerequisites 

Before considering the rotational and stretching properties of bilipschitz 
and quasiconformal mappings, we briefly recall the basic concepts underlying 
our study. 

2.1. Bilipschitz and Quasiconformal Maps. By definition, in any di- 
mension n ^ 2, X-quasiconformal mappings are orientation preserving home- 
omorphisms / : O — > Q! between domains £1,0,' C M. n , contained in the 
Sobolev class W^(Q,), for which the differential matrix and its determinant 
are coupled in the distortion inequality, 

\Df(x)\ n K det Df(x) , where \Df(x)\ = max \Df(x)£\, (2.1) 
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for some K ^ 1. Prom now on, we will consider these mappings only in 
dimension n = 2. It is clear that any orientation preserving L-bilipschitz 
mapping is ivT-quasiconformal with K = L 2 . 

In comparison, any quasiconformal mapping of the entire plane is qua- 
sisymmetric, i.e. satisfies the estimate 

\f(x)-f(z)\ (\x-z\\ „ , , 
, x,y,zeC, (2.2) 



\f(y)-f(z)\ ^ '\\y 
for some continuous strictly increasing n : M+ — > R+ with rj(0) = 0. We will 
make frequent use of this geometric description. Conversely, any quasisym- 
metric mapping is iT-quasiconformal with K = 

What makes quasiconformal mappings particularly flexible in dimension 
n = 2 is that they satisfy the Beltrami equation 

K - 1 

fz = Vfz with ^ k < 1, k = (2.3) 

A + 1 



1 2 \ k 

Conversely, any homeomorphic W,' -solution to (2.3) is A'-quasiconformal 



with K = (1 + k)/(l — k). We refer to [2] for definitions and basic facts on 
planar quasiconformal maps. 

If the coefficient [i is compactly supported, the mapping / is conformal 
near oo and we may use normalisation f(z) = z + o(l) as z — > oo. In 



this case we call / the principal solution to (2.3); such a solution is uniquely 
determined by the Beltrami coefficient fj,(z) and it is a homeomorphism. The 
Beltrami equation paves the way for embedding the principal solution / into 
a holomorphic family of quasiconformal maps. One may simply consider the 
flow of principal solutions {f X (z)}\<=B to the equation 

h = H\fz, »\(z) = A/i(z)/|M|oo, f = f X - (2-4) 



Similarly, the global solutions to (2.4) normalized by /(0) = 0, f(l) = 1 de- 



pend holomorphically on A G D. More generally, within both normalizations 
we achieve the holomorphic dependence as soon as ^\{z) varies holomorphi- 
cally with the parameter A. For details and further facts on holomorphic 
dependence see [H Section 5.7]. 



2.2. Choosing the logarithmic branches. In describing the different as- 
pects of rotation under quasiconformal mappings we will need a convenient 
and systematic way to discuss the various branches of the logarithm of the 
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difference f(z) — f(w). The following simple observation serves best our 
purposes. 

Proposition 2.1. Let f : Q — > M? be a homeomorphism of a simply con- 
nected domain ficR 2 . Assume that f is differentiable at a point wq G Cl 
with positive Jacobian. Then the logarithmic expression 

log - , (^)e!!x!l, z^w, (2.5) 
z — w 

admits a single-valued continuous branch. 

Moreover, any continuous branch is uniquely determined by its value at 
any given pair (z,w) G SI x Q, z ^ w. 

Proof. Since / : f2 — > M? is differentiable at wo G and J(wo,f) = 
\fz(w )\ 2 - \fz(w )\ 2 > 0, the function 

z f{2) ~ f{Wo) (2.6) 
z - w 

is continuous and non- vanishing in fi\ {wq}. Its accumulation set at z = wo 
is contained in a closed disk B(a,r), centered at a = f z (wo) and of radius 
r = \fz(u>o)\ ^ . By our assumptions on the Jacobian, £ B(a, R) as soon 
as R > r is sufficiently close to r . If 7 is any loop in 0\{it?o}, it is homotopic 
to a loop that lies arbitrarily close to wq, i.e. to a loop whose image under 



/ lies in B(a,R). Therefore (2.6) admits a single valued logarithm as a 
function of z 6 £1. 

In order to show that log ^ z ">~^ w ^ admits a single- valued continuous 
branch in the whole domain U := (f2xS2)\{(w;, w) : w G f2}, it is then enough 
to verify that any loop in U is homotopic to one that lies in the section 
U n {(z, w) : w = wo}. Since is simply connected, it is homeomorphic to 
the unit disc D and we just need to consider the case U = D . Moreover, we 
may assume that wo = 0. Given any loop (a, j3) : [0, 1] — > D X B that avoids 
the diagonal (that is, a(t) 7^ (3(t) for allO ^ t ^ 1 ), the required homotopy 
is given by 

Indeed, a s {t) / f3 s (t) for all < t < 1 and < s < 1 . □ 
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Remark 2.2. In most cases that we will encounter there is a natural choice for 



the branch of the logarithm (2.5 ). For instance, if / : C — > C is normalized by 
/(0) = 0, /(l) = 1 we will choose log/(l) = log ((/(l) - /(0)/(l - 0)) = 0. 
Or if / : C — > C is a principal solution to (|2.3[) with f(z) = z + o(l) as 



z — > oo, then we consider the continuous branch with 

f(z) - f(w) 

log — » as z — y oo. (2.7) 

z — w 

and call it the principal branch. 

Remark 2.3. Later on, we also need to consider flows of homeomorphisms 
f x : C — > C that depend continuously on a complex parameter A 6 D. 



Proposition (2.1) generalizes to this situation, and there is a single-valued 



continuous branch of the logarithm 
f x (z) - f x (w) 

log — v ; defined for (z,w) € C x C, z^ w, and A € D. 

z — w 

In order to prove this one observes that any loop 7 = (71,72,73) in the 
domain 

((C x C)\{(z,w) : z = w}) x D 

is homotopic to a curve lying in ((C x C) \ {(z,w) : z = w}) x {0} by the 
trivial homotopy (t, s) h-> (71 (t), 72(^)1 (1 — s )73(*))- This reduces one to the 



case that was already handled in the proof of Proposition (2.1). Naturally 



the branch is uniquely determined by its value at any (z, w, A) with z ^ w. 



3. Notions of rotation 

We start with the pointwise notions, i.e. describe the extremal behavior 
and the optimal bounds for the rotation and stretching at a point zq £ C. 

There are (at least) three natural different and geometric ways to describe 
the rotational properties under a planar mapping: the infinitesimal, the local 
and the global concepts. They have different and complementary geometric 
and analytic descriptions, but we shall see later that they are all intimately 
related. These relations will then form the basis of the quasiconformal and 
bilipschitz multifractal properties studied in Sections [5] and [6| 

We begin with the local point of view, where one fixes the argument at 
a chosen pair of points zq 7^ z\ G C, and then studies the behavior of the 
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logarithmic cross ratio 

m - f(z ) 



log 



as z 



z , 



Z / Zq. 



J(zi)-f(z ) 

The scale invariance of this expression allows universal distortion bounds. 

As a next step, the local bounds enable one to study the geometric rate of 
spiraling for the image of an infinitesimal line segment; for precise definitions 



see (3.8) - (3.10) and Theorem 3.3 



Finally, for principal mappings; that is, when f(z) = z + o(l) as z — > oo, 
we have the third possibility of taking z = oo as the reference point. Here it 
is natural to travel from oo to zq along a line segment and study how much 
the image rotates around f(zo). 

3.1. Rotation and Stretching: The local picture. By the classical 
Mori's theorem, a iT-quasiconformal mapping / is locally Holder contin- 
uous with exponent 1/K. We look for a similar bound on the maximal 
change in argument. It turns out that a fruitful point of view is to measure 
twisting and stretching simultaneously, and this is one of the leading insights 
in what follows. By scale invariance it makes no difference to take zq = 
and z\ = 1 and assume f(0) = 0, /(l) = 1. 



Theorem 3.1. Consider a K -quasiconformal mapping f : C 
ized by /(0) = and /(l) = 1. Then for any < r < 1, 



normal- 



log f(r) 



K + 



K 



log r 



1 



K ~jA l °S 1 + <K), (3.1) 



where the constant c(K) is independent of r. Here the continuous branch 
of the logarithm log/(r) is determined by log/(l) = 0. 

Furthermore, for any < r < 1 and for any number r on the circle 



5 M 



(3.2) 



the equality in (3.1) holds with c{K) = for the K -quasiconformal mapping 
/rW = nN T , r = a(l + i 7 ). (3.3) 



For an illustration of the parameter disk (3.2), see Figure [2] presented 
after Theorem 13.31 



10 



K. ASTALA, T. IWANIEC, I. PRAUSE, AND E. SAKSMAN 



Before embarking to the proof, let us recall some of the earlier results on 
quasiconformal rotation. Gutlyanskh and Martio [T2] answer a problem 
of John [15] on bilipschitz mappings by determining the maximal rotation 
for K-quasiconformal maps that fix given annuli, keep one boundary circle 
fixed and rotate the other boundary circle. The extremal in this situation 
is given by the map (3.3) where a = 1, corresponding to a pure rotation. 



Balogh, Fassler and Platis [I] extend this result to annuli with different 



modulus, and the extremal is of the general form (3.3). Both of these works 



actually consider a fairly general class of maps of finite distortion. However, 
they only consider mappings between round annuli. Our result relaxes on 
this hypothesis, and we simply ask for the maximal combined rotation and 
Holder distortion of f(z + r) — f(z) for any mapping defined in the entire 



plane; after a normalization we are reduced to (3.1). 



Proof of Theorem 3.1, The case of the equality for f T is immediate. Note, 



in particular, that the circle (|3.2|) is precisely the set of points r such that 

K - 1 
K + 1 ' 

f T in (3.3) satisfies the Beltrami equation fz = ^j|/ z , the 



r - 1 
r + 1 



Since each / 



mappings are indeed i^T-quasiconformal. 

Next, let < r ^ 1, and consider first iT-quasiconformal mappings which 
are conformal outside the disk B(0, 2) with 



f(z) = z + 0(l/z) as \z\ — > oo. 



(3.4) 



Embed then / to a holomorphic flow {/a}asd as in (2.4), so that each map 



has the same asymptotics (3.4). By Koebe distortion (see e.g. [21 Theorem 



2.10.4]) f\B(0, 2) C B(0, 4), so that the values of f\(z)-f\(Q) for \z\ < 2 all 
stay within the disk B(0, 8). We thus observe that the holomorphic function 



V>(A) := log 



A € 



/a(0-/a(0)' 

where the branch for the logarithm is fixed by setting "0(0) > 0, has positive 
real part. 

These preparations take us to the main point of the argument, where we 
simply apply Schwarz lemma to ip : D — > H, with H the right half plane. 



This implies 



1+|A| 2 

HAp 



2l Al 

^(0) ^ i-|A| 2 ^(0)- Unwinding the definitions 
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nves us 



log[/(r) - /(0)] 



K + — ) log r 



(3.5) 



Thus a version of (3.1) is shown to hold for the principal branch. 



It remains to find estimates when we do not have conformality outside 
2D. For this, given a quasiconformal mapping / of C normalized by /(0) = 
and /(l) = 1, apply Stoilow factorization ([2 Theorem 5.5.1]) and represent 
/ = F o / , where /o is a principal map conformal outside the disk 2D, and 
F a quasiconformal map which is conformal in the domain f2 = /o(2D). 

Here we will make use of the following lemma. 

Lemma 3.2. Suppose that F: f2 — > C is conformal, where $7 C C is simply 
connected. Then for any x,y,w G Q with x ^ w, y ^ w, and any chosen 



branch of the logarithm {c.f. Proposition 2.1) it holds that 

F(w) . F(y)-F(w) 



log 



Fx) 



log 



^ 10pn(x,y). 



x — w y — w 

where is the hyperbolic metric of the domain Q. 

Proof. The lemma is basically a reformulation of the classical distortion 
bounds for conformal mappings in the unit disk. Indeed, suppose first that 



F is conformal in D with -F(O) = 0. Then |19l p. 21] we have 



6(1- 
that 



1 uniformly in D which implies for any branch of the 



F"(z) 



F'{z) ; 
ogarithm 



\\og F'(x)- log F'{y)\ ^3 Po (x,y), 
On the other hand |19| p. 66], we also have 



x,ye 



log 



zF'(z) 



«S log 



1 + 1*1 



z G 



F(z) "~~°l-\z\ 
with the normalization log 1 = at z = 0. Combining the estimates gives, 
again for any branch of the logarithm, 



, F(x) , F(y) 
log - log 



x y 

In the general case compose with the Riemann map 
Then for any choices of the respective branches 



x,y G 



(3.6) 

n, 0(0) = w. 
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f(0(z))-f(0(o)) Fj^z)) - nm) , Hz)-m 

8 4>{z) -0(0) 8 z 8 z 

+ 2vrin, 



for some n£Z. Applying now (3.6) twice proves the claim. □ 



Returning to the proof of Theorem |3.1| and the factorization / = F o / , 
the triple x = fo(r),y = /o(l), w = fo(0) is contained in fo(JS>), and this set 
has hyperbolic diameter in f2 = Jq(2D) bounded uniformly in terms of K 
only. As F(x) = f(r),F(y) = 1 and F(w) = 0, Lemma |3 . 2| shows that 



log/(r) -log(/ (r) - / (0)) +log(/ (l) - / (0)) ^ M(K) < oo. 



On the other hand, applying (3.5) with r = 1 to /o implies 



log(/ (l)-/o(0))|< (if -1) log 8. (3.7) 



Finally, combining these bounds with (3.5 ) for a general < r < 1 completes 



the proof of the theorem. □ 

3.2. The Infinitesimal Rate of Rotation. For a quasiconformal mapping 
/ : fi — >• fi' defined in a proper subdomain C C it is senseless to ask for 
uniform bounds for rotation or stretching in all of £1 - the general properties 
are determined by the geometry of the specific domains in question. One 
should recall in this context that even a conformal map from unit disc onto 
e.g. a snowflake domain rotates in an unbounded manner as one approaches 
suitably the boundary. 

On the other hand, in any domain for infinitesimal rotation and stretching 
we have uniform bounds. To describe this more precisely, first consider for 
t := a(l + ry) the prototype example 

/ r (z) = A| z |«(i+n) 5 for |z| ^ 1 with f T (z) = z, for \z\ > 1. (3.8) 



The map f T takes a radial segment in the unit disk to a logarithmic spiral 
with rate of rotation 7 G M, see Figure [TJ while a > can be considered as 
the stretching exponent of f T at the origin. 
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Figure 1. A logarithmic spiral with rotation rate 7 = 5. 

Intuitively, keeping the model map f T in mind, for a general mapping / 
the local exponent for stretching at a point zq G should then be given by 

log 1/(^0 + 0-/(^)1 

«(«>) = ^ ' ( 3 - 9 ) 

while the (geometric) rate of rotation one would define as 

/ \ r arg(/(z + r) - f(z )) 

7(2:0) = hm — — — -p. (3.10) 

r->o log |/(2o + r)-/(2:o) I 

In defining the latter expression one takes into account that the more the 
mapping is compressing the more it has to rotate to produce a given geo- 
metric spiral as the image of a half-line emanating from zq. An easy way to 
illustrate this is to consider the curve t — > t 1+t "< and the change of variable 
t = s K . 

The problem with the above intuitive notions is, of course, that the re- 
spective limits need not exist as / may have different behaviour at different 
scales. Therefore we formulate the infinitesimal joint stretching-rotation 
bounds as follows. 

Theorem 3.3. Assume f : $7 — >■ f2' is a K -quasiconformal homeomorphism 
between two planar domains, and that zq € £1 is given. 
Suppose we can find radii r^ — > 0, such that the limits 

a = hm 1Og ' /(Z ° + rfc) - /(Z0) ' , (3.11) 

fc-S>oo log Tfc 

,. arg(/(2 +r fc ) - f(zp)) 

7 = nm -j — 777 , \ Ti 77 {3.12) 

fc^oc log /(20 + r k ) - f(zo)\ 
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Remark 3.4. To explicate the statement above, in ( 3.12[ ) we use some con- 
tinuous branch of argument r i— >■ arg(/(zo + r) — f(zo))- However, it is 
important to notice that the limit ( |3.12[ ), if exists, is independent of the 
choice of branch used. 

Corollary 3.5. In particular, for every K-quasiconformal mapping 
f : O — > O' and for every point zq £ each limiting rate of rotation 7 



in (3.12) satisfies 



|7K 



1 



K 



1 

K 



(3.14) 



The tangent line in Figure [2] illustrates the extremal case in Corollary 3.5 

= the model map (3.8) is K- 



Recall furthermore that for 
quasiconformal, that for zq = the exponents a and 7 are given by (3.11) 



and (3.12), respectively, for any sequence of radii — > 0, and that the 



exponent r = a(l + £7) satisfies (3.13) as an equality. Thus the bounds of 
Theorem |3.3| are optimal. 



Proof of Theorem 3.3. It suffices to consider simply connected domains 
ficC. We first apply Stoilow factorization [2, Theorem 5.5.1] and represent 

f(z) = ipof (z), zeQ, 



QUASICONFORMAL MULTIFRACTAL SPECTRA 



15 



where /o : C — > C is K— quasiconformal and <p is conformal in the domain 
n" = / (n). If = C we take </?(z) 
Theorem |3 . 1 1 now implies 



z. 



log 



fojzp + tr ) - fo(z ) 
fo(zo + r ) - f {z ) 



(3.15) 

te (0, l], 



where ro is such that B(zo,ro) C Vt. On the other hand, Lemma 3.2 gives 



log- 



log 



y 



^ 10pn//(x,y). 



for x,y,z € O". Choosing a; = fo(zo + iro), y = fo( z o + r o) an d -2 



fo(zo) shows that (3.15) holds for /, too, with c(K) replaced by a constant 



depending on K and the hyperbolic distance in Q" between fo(zo + tro) 
and fo(zQ + ro). Adjusting ro this distance can be made arbitrarily small 
uniformly in t € (0, 1]. 

Finally, letting t = r^/ro, dividing by log(l/rfe) and taking k — > oo gives 



the bound (3.13). Note here that by Mori's theorem we necessarily have 



3.3. Global Aspects of Rotation: The Complex Logarithm log/ z . 
With Theorems |3.1| and |3.3| one has a first description of the rotation and 



stretching properties of quasiconformal mappings. Theorem 3.3 determines 
the extremal infinitesimal phenomena, and it is of interest only at the points 
of non-differentiability 

For a principal mapping with f(z) = z + o(l) when z — > oo, it is also 
natural to consider a global notion, the rotation the image of a half line 
from oo to zq makes around the image point /(zq), see Figure [3j This global 
measure of rotation, however, works best at points of differentiability and 
as it turns out, is described by the function log/ z . 

Before embarking into the several interesting analytic and geometric prop- 
erties that the complex logarithm of f z encodes, we of course need to clarify 
the pointwise definition of the function log f z . There are actually two al- 
ternative ways for this: an analytic approach and a geometric one. 
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Analytic definition. As explained in Subsection 2.1 above, by letting fi(z) = 
H\(z) depend holomorphically on a parameter A £ D, we have natural ways 
to embed a given principal mapping / into a family of quasiconformal maps 
f x , with analytic dependence on the parameter A £ D. To make further use 
of this property, we need a refinement of [31 Lemma 3.7]. 

Lemma 3.6. There is a set E C C of full measure, such that for every 
A 6 D and for every z E E, the map / A is differentiate at z with f£ ^ 0. 

In addition, for a fixed z £ E and for any compact icD the differential 
quotients 

/A(z + t ;- /A(z> . t eM ,XeA, 

are uniformly bounded. 

Furthermore, for z G E the function A \-¥ fz(z) is holomorphic in A £ D 



and f/ie equation (2.3) ura'i/i fJ. = fJ.\ holds true pointwise in E. 



We postpone the proof of the Lemma to the end of this section. As f°(z) = z, 
it is natural to set logf z = 0, and thus requiring that log f z (z) remains 
holomorphic in A defines the logarithm uniquely at every z £ E. Returning 
to the original solution / we can now set the analytic definition 

log f z (z) = log f^(z), k = Halloo, ze E. 

Geometric definition. While the analytic approach to log f z requires map- 
pings defined in the entire plane, there is a more geometric approach that 



works in general settings. Indeed, by Proposition 2.1 given a quasiconformal 
mapping / : Q — > Q' between simply connected domains, each branch of the 
logarithm 

log *W ~ /H ( ZjW ) € fixn, z^w, (3.16) 
z — w 

determines at points of differentiability of / (such that also fs(z) = li{z)f z (z), 
hence almost everywhere) a corresponding branch of the complex logarithm 
of f z , namely 

■f(z + t)-f(z) 



log(9/(z)) = lirn log 



log[l + /i(z)]. (3.17) 



t 

The role of the auxiliary horizontal approach is of course redundant, and 
one could as well apply an approach from any other direction, 
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log (df{zj) = liinlog 



log[l + e- 2l V(^)] 



For general / : O — >■ O' there is no natural choice for the argument of 
df(z) in (3.17) but for principal mappings / : C — > C, unless otherwise 



explicitly stated, we will always consider the principal branch of log f z (zo), 



determined by the conditions (2.7) and (3.17). This branch admits a clear 



geometric interpretation. Namely, we may consider the corresponding con- 
tinuous branch of the argument, 

-f( Zo + t )-f(z y 



arg[/(z + t) - f(z )] = Imlog 



t 



< t < oo, 



where in the case of the principal branch ( |2.7[ ) we have arg[/(zo + t) — 
f( z o)] — > as t — > +oo. The integer part of ^ arg[/(2:o + to) — f( z o)] then 
gives the number of times the point f(z) = f(zo + t) winds around /(zq), 
when z = zq + t moves from t = +oo to t = to along the horizontal line 
through zq. Further, at any point zq of differentiability of / the argument 
satisfies 

-f( Zo + t )-f(z y 



axg[f(zo+t)-f(zo)] 



arg 



& Tg[df(z )+K z o)df(z ) + 0(t)] 



as t — > + . Since ^ k < 1, we may choose arg[l + [a(z)] £ (— f,f), 
and obtain almost everywhere the geometric definition of the argument 
arg[<9/0o)] = Imlog df(z ). 

As the above discussion indicates, with this definition at points of dif- 
ferentiability, for the principal branch we have | Im \og(d f)(z) — 2irn\ < tt 
where n is the number of times f{w) winds around f(z) when w travels from 
oo to z along a fixed radius, see Figure [?| 

Even if above we used two quite different approaches to the complex 
logarithm of f z , both methods lead to the same concept. 

Lemma 3.7. For any principal quasiconformal mapping f : C — > C, the 
geometric and analytic definition of the principal branch of log f z agree on 
a set of full measure. 
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Figure 3. argdf(zo) measures the total winding of the im- 
age curve around the point f(zo). 

Proof. It is useful to observe that the definition of the principal branch of 
arg[/(zo + t) — f(zo)~\ given above can equivalently be obtained from the 
imaginary part of the analytic continuation /i(A) = log(/ A (zo + 1) — f x (zo)), 
when one sets h(0) = logt. 

For the proof we only need to consider the set £cC given by Lemma 



log[l + /Z A (*)] 



3,6 At any z G E define 

ot(A, z) = log 

Then for a fixed z G E, limt_>.o+ (h(X, z) is holomorphic in A, as a pointwise 
limit of bounded holomorphic functions (here we use the statement of uni- 
form boundedness of the differential quotients in Lemma |3.6[ ). Further, the 
function vanishes at A = and one has exp (lirm : _ > o+ ctt(A, zyj = df x (z) for 
every A G D. □ 

In particular, the definition of log f z does not depend on how the mapping 
is embedded to a holomorphic flow. 

Once the pointwise notion of log f z is settled, one can ask for its function 
space properties. Here it follows from the work of D. Hamilton [13] that 
log f z G BMO(C). He complexifies an argument of Reimann |21j and esti- 
mates the A-derivative of the function. Indeed, writing f x+€ = h e o f x and 
using [2, Thm. 5.5.6] we have 

N^wfqwf + °( e2 )= : ^°/ A (*) + 0(e>) 

(3.18) 

For the derivatives of h e we have at almost every point (h e ) z = l + eS(v\) + 
C(e 2 ) with (h e )z = tv\ + C(e 2 ), see [2J Section 5.7], so that from the chain 
rule 
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d z f X+e - d z f x = d z h e (f X ) d z f X - d z f x + d- z h e {f x ) d- z f* 



= e(Sv x )(f x )d z f x + ev x (f x )d z f^ + 0(e 2 ) 

Thus 



° yr} -{Svx){f X ) + vx{f X )lTx-4r = (^a)(/ A ) + A 



Since by another theorem of Reimann, quasiconformal mappings preserve 
the space BMO(C), the right hand side has £?MO-norm uniformly bounded 
for |A| ^ fco < 1. By integrating along a radius from A = 0toA = || / u|| oo we 
obtain 

Proposition 3.8. |13] Suppose f G W^i^) is a principal quasiconformal 
mapping. Then 

logf z EBMO(C). 

In particular, arg f z is exponentially integrable, i.e. 

e 6 ' arg ^ ' E L\ oc for some positive constant b. 
One of the key points of the present work is finding the optimal form of the 



exponential integrability. This will be done in Corollary 4.10 where we give 
the sharp bounds for b in terms of the distortion K(f). 

To complete the subsection we need to justify Lemma (3^6} 



Proof of Lemma 3.6 After [3, Lemma 3.7] one only needs to verify the (si- 
multaneous and uniform for all A) differentiability almost everywhere. Recall 
that in plane continuous functions / £ W^(C) with p > 2 are differentiable 
almost everywhere. A proof is based on the classical Morrey estimate (see 
[8, Section 4.5.3]), 

\9(y)-9(x)\ < c p \y-x\ ( 1 f \D 9 A , g 6 W%(C), 

ytS{X, \y — X\) JB(x,\y-x\) J 

where p > 2. The differentiability a.e. follows by applying the above estimate 
to the function g(y) = f(y) — f(x) — Df(x)(y — x) and using a Lebesgue 
point argument. 
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For our purposes we only need to quantify this a little bit. Denote by 
Qj(x) the maximal difference quotient at x, given by 

« /W := sup mjs. 

ly-xjsgi \y ~ x \ 

As a direct consequence of Morrey's estimate, any continuous / G W,J!!(C) 
satisfies the pointwise estimate 

Q f (x) < c p M p (Df)(x), 

B^F) /s(x,r) lfl'l P ) • Es P eciall y> wnen deriva- 
tives of / belong to L P (C) the weak (l,l)-continuity of the maximal function 
yields for any s > 

|{Q />S }K4 S ~ P P/U£- (3.19) 

In our situation, for any given 5 < 1 we may pick p = p{5) > 2 so close 
to 2 that the standard ££ v -properties of the Beurling operator and the 
Neumann series representation (see [21 p. 163]) enable one to write, for any 
A G D, 

oo 

f x (z) -z = J2 ^fn(z) with [|D/ n || p < co 5- n ' 2 . (3.20) 
n=l 



Here cq depends only on the size of the support of \i in (|2.4|), and each 



fn G C(C)n^ f(C) for every q < oo. The series converges locally uniformly 



m <u 

For n ^ 1 write := ££L n A fe /fc • Fix <5 < 1, denote A k := {Q fk > 
<r 3fc / 4 }, and use ( |3.19[ ) to estimate |A fc | < eg 5-^/2^/4 = C( 5fc P /4_ Next> 
write F n := \J™ =n A k , so that |F n | < 5 np ^ while if a? G F£ we have 

00 

Qva(*) < E 5fc(5 " 3fc/4 £ 571/4 for ever y i a i < 5 - ( 3 - 21 ) 

fc=n 

Clearly this implies that the set E$ := n^ =1 F n has full measure with 
lim^oo Q hnX (x) = for every |A| < 6 and x G F 5 . 

Pick now a set £7 of full measure such that each f k in (3.20) is differen- 
tiable at every point of E. Especially, for every |A| ^ S the function f x is 
differentiable at each x G E n Eg, with 



Df\x) = Id + J2^ k Dfk(x). 



k=l 
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Whence we may choose E := EnU^j ^i-i/e- Finally, the statement about 
uniform boundedness of the differential quotients with respect to A in com- 

□ 



pacts follows directly from (3.19) and (3.21). 



4. Interpolation with complex exponents and Burkholder 

integrals 

4.1. Interpolation with complex exponents. Let (O, a) be a measure 
space and let ^(0, a) denote the space of complex-valued a -measurable 
functions on Q. We consider Jz? p ($l,a) spaces in which the (quasi-)norms 
are defined by 



1*1 



< p < oo , and ||<J>| 



ess sup |$(z) 



We shall consider analytic families &\ of measurable functions in O, i.e. 
jointly measurable functions (x, A) i— >■ Q\(x) defined on f2 x U, where U C C 
is a domain, and such that for each fixed x G the map A i— > <&\(x) is 
analytic in U. The family is said to be non-vanishing if there exists a set 
E C 0, of a- measure zero such that &\(x) ^ for all x G f2 \ E and for all 
AG U. 

In the following we give a generalization of the interpolation results from 



[3j to complex exponents. In fact, restricting below in (4.1) to the real 
interval /3 G po, ^ f 1 + p^] one arrives at the bounds of [31 Lemma 1.6]. 



Lemma 4.1 (Interpolation Lemma with complex exponents). Suppose 
{&\ ; |A| < 1} C cr) is a non-vanishing analytic family of functions, 

parametrized by complex numbers A £ D, such that for some po > 0, 

= 1) a nd II^aIIpo ^ 1 f or ev ^ r D A G B. 
Then, for every |A| < 1 and for every complex exponent (3 G C contained 
in the closed ellipse 



+ 1/3 -Pol < 



Po_ 
|A|' 



we have 



J 

Jn 



do- sC 1. 



(4.1) 



(4.2) 



The choice of branch in (4.2) is the natural one, determined by the condition 
log<I>o = 0. 
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Proof. By considering the analytic family cj>^ 0//2 we may restrict our atten- 
tion to the po = 2 case, which will later make calculations simpler. Observe 
that our assumption ||<&o||2 ^ 1 implies <r(f2) ^ 1. Actually, by the maxi- 
mum principle for analytic Jf 2 - valued functions, we may further assume the 
strict inequality 

<r(fi) < 1. (4.3) 

Otherwise <3?a would be constant in A, as we have 1 1 a 1 1 2 ^ 1 for ah A £ D. 
We may also assume in the proof that < c ^ | < 1 ) a( 2; )| ^ C < oo uniformly 
for all (x, A) G f2 x U, as the reduction of the general situation to this is 
done exactly as in [31 Section 2]. Similarly, we choose an arbitrary positive 
probability density p, uniformly bounded away from and oo, 



G^#(H, cr), ||p ||i = / p(x)dx 



By Jensen's inequality using the convexity of x i— > xlog(x) and (4.3) we 
have I := f n p(x) log p(x) dx > 0. Temporarily assuming that p is fixed, 
we consider the holomorphic function / in the unit disk 

fW = 7 / p(x)log$ A (x)dx . 



I . 

Again by Jensen's inequality we have the bound 

2Re/(A)-l = j [ p(x)log^4^ d ^ 7 lo g(7 ^ ^TV" dx ) 
I Jn P(x) I \Jn P{x) J 

= j log||$ A ||| <0, for |A| < 1. 

Thus / maps the unit disk into a half-plane Re /(A) 1/2, while /(0) = 
by our assumption = 1. At this stage we appeal to Schwarz lemma and 
deduce that for any ^ k < 1, the image /({|A| ^ k}) lies in a hyperbolic 
disk Dk, centered at 0, of the above half-plane. Precisely, 

f({\\\^k})cD k = {-^:\z\^k 



J + z 

Our objective is to find all exponents j3 G C such that for all |A| $C k, 

Re ^ A )) = TTTT Re (P [ Pip) log $a0*0 dx) < 1. (4.4) 

Equivalently, we aim at the estimate 
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and once this is achieved, by specialising the choice of p, that is, choosing 

-l 



p(x) := ISaWI 



, we will obtain our assertion 



log 



da ^ 0. 



Given |A| k, the condition ( |4,4[ ) is now ensured by 

Re(f3w) ^ 1, for all w 6 Djt, (4.5) 
as the range of /(A) lies in D^. To state this requirement more explicitly, 



note first that as 



z 

l + z 



l + z 



l + z 1 



inequality (4.5) takes the form 
l + z' 



Re[P z 



l + z 



z\ 2 Re/3 < 1 



\z\ < k. 



(4.6) 



Here only the first term depends on the argument of z, with z 
preserving the circle of radius k. Thus (4.5) is equivalent to 

k \P\ - k 2 Re < 1 - k 2 

By simple algebra, we have here the equality if and only if 

1 J\ — £2 
P = l+ cos9 + i^—. sine, for some 6 G [0, 2tt] . (4.7) 



Thus the extremal P lie on the ellipse with foci {0,2} and eccentricity k, 
so that (4.5) is equivalent to k\P\ + k\P — 2\ ^2. For a general exponent 
Po > 0, the ellipse takes the form of (4.1). □ 



Remark 4.2. Given P, the left hand side of (4.6) attains its maximum over 
{A : |A| = k} at the point where 



that is when 



P 



arg/3 = 2 arg(l + A) — arg A, 
1 + A 



A 



(l + A)s(A), s(A) € M-| 



Testing this requirement against ( |4.5[ ) shows that s(A) ^ 1/(1 + Re A). 
Therefore the equality in (4.5) is attained at A, |A| = k, if and only if 



P 



(1 + A) 2 
A(l + Re A) 



1 + 



1 



1 + A 



A / 1 + ReA 



(4.8) 
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4.2. Burkholder integrals. Let /: C — > C be the principal solution of a 
Beltrami equation 



f- z (z) = p(z)f z (z), \p(z)\^k X »(z), 0O<l- 



(4.9) 



In studying the rotation spectrum of bilipschitz or quasiconformal mappings 
we are faced with the question for which exponents /3 G C is the complex 
power (f z ) locally integrable? It turns out that the universal bounds are 
given exactly by (4.1) with p = 2, |A| = k, that is, in terms of an ellipse 
having foci {0, 2} and eccentricity determined by the ellipticity constant k 



of the equation. See Theorem 4.9 for the precise statement. 

In fact, we are going to carry out our analysis in the weighted setting and 
consider the so called Burkholder type integrals. Here recall the functionals 
introduced by Burkholder [7], which applied to the derivatives of a map 
/ : M 2 ->■ M 2 take the form 



BJDf) 



^(pJ(zJ) + (2-p)\Df\ 2 



\fz 



■ \D/r 2 

( P -i)\M)-(\f z \ + im r 1 , 



Originally the functional was discovered by Burkholder in his studies of 
optimal martingale inequalities - since then optimal integral identities re- 
lated to Bp, and in particular its conjectured quasiconcavity (for |jp — 1| ^ 1) 
have been of wide interest. For recent advances and background see e.g. [3], 

Here we search for corresponding functionals determined by a complex 
parameter p G C \ 5(1,1). We define an auxiliary unimodular function 
p = p(z), \p\ = 1, by requiring that the complex numbers 



pp(z) and 1 + p(z)\p(z)\ have the same argument. 



(4.10) 



Theorem 4.3. Suppose we are given a complex parameter p with 1 ^ \p ■ 
1| ^ J and an exponent (3 G C with 



+ |/3-2| <2|p-l|. 



Then for every principal solution to (|4.9|) we have 
1 



7T JO 



\M+P\fz\ ~\p\\fz 



fz+ P\p\fz 



< 1, 



(4.11) 



(4.12) 



where the unimodular function p = p(z) is determined by (4.10). 
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Figure 4. The relation between p and /3 in (4.13), 
\[3\ + \p - 2| = 2\p - 1| and Re(/3/p) = 1. 

The estimate holds as an equality for f(z) = z. Furthermore, when f3 is 
determined (uniquely) in terms of p by the equations 

| J 9| + | / 9-2|=2|p-l| and Re(/3/p) = 1, (4.13) 



we have the equality in (|4.12|) for every power map of the form 

z 



77 

where p G [0, 1] and Id = k. 

1 + r? 



In other words, the functional (4.12) is quasiconcave for the parameter 
values 1 ^ \p — 1| ^ 1/fc, within the class of principal /c-quasiconformal 
deformations. 

Remark 4.4. Using the disk filling procedure as in [3], we get many more 
extremals for the functional (4.12). 

Remark 4.5. For p 2 we have p = 1, while for p < 0, p = —1. Thus for real 



p with the choice of (4.13) we get back the Burkholder functionals B p {Df). 
On the "phase transition" boundary \p— 1| = 1, (4.13) forces (3 = 2 and we 
recover the Jacobian. 



Proof of the Theorem J^3_. The proof adopts ideas from [3] to the case of 

complex exponents. Given the Beltrami equation (4.9), we define a holo- 

morphic variation as follows. For A 6 D, set 

— - p(z) a x {z) p(z)\p(z)\ A 

p\{z) = ax{z) p(z) ■ -—— where — —=p-— — — — - • — — 

\p(z)\ l + a\(z) 1 + p(z)\p(z)\ 1 + A 

Above we use the convention "0/0 = 0" whenever dividing by zero. 
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Our assumptions ( 4.10| ) and \p — 1| ^ | show that p ■ j^jjj G [0,1] 
and thus | /^a C- 2 ) I = | c^a (^) I ^ |A| < 1. This makes legitimate to solve the 
Beltrami equations 



under the normalization of the principal solution. 

We recover the original equation for the complex value A = and the 
Cauchy-Riemann equations for A = 0. 

Next, interpolate the analytic family of functions given by 

d> A (z) = (l + a A (z))/, A (z)/0. 



Indeed, <l?o = 1 an d according to Lemma 3.6 {<3?a} is a non- vanishing family 
in the sense required by Lemma |4.1| Furthermore, we have the comparison 



1 - 2 Re 



1 + OL X {z) 



p- 



l + p(z)\p(z)\ 



--: u(z). 



Thus by the classical area theorem, see e.g. [21 p. 41], we have the Jz? 2 -bounds 

' 1 '*a| 2 w< 1, 

with the measure space ^#(D, dz) to 



7T 



and we may apply Lemma 
obtain 



4.1 



1 



P 



fz+p\p\fz 



< 1, 



(4.14) 



with G C as in (4.11). Since by (4.10) the complex numbers p p(z)\p(z)\ 
and 1 + p{z)\p{z)\ have the same argument, the integrand, in fact, takes the 
equivalent form of (4.12). 



Concerning sharpness, let f(z) 



\z \ 1+r i with 



k. As pf(z) 



IV, 



Since f z 



(4.12). 



^ determines the unimodular factor p = r]/\r]\, 
1 one computes that Re = Pf^- 
jq— |z| , a direct substitution shows that the equality holds in 

□ 



the requirement pj 
i.e. r\ = p\p\. From Re(/3/p) 

-877 



4.3. Higher complex integrability. By Stoilow factorization, for local 
integrability issues it is enough to control the behaviour of principal maps. 
Hence, our previous theorem yields immediate corollaries. 
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Theorem 4.6. Suppose f : C — > C is a K -quasiconformal mapping and 
B = B(z,r) C C is a disk. Then for any exponent (3 G C such that 

K + l 



+ |/3-2| < 2 



K-V 



we have 



f(z + r)-f(z) 
r 



1 

W\ 



B 



ft 



(4.15) 



(4.16) 



f(z + r)-f(z) 
r 



where the constants c\,02 depend only on K and j3 



Remark 4.7. The branches of log f z {w), for w G B, and of 
are here chosen as in (3.16) - (3.17), using the same branch for both. More 
precisely, consider any of the branches of log ^M^^i discovered in Propo- 



sition 2.1 As explained in (3.17), this determines at a.e zeCa branch of 



log df(z), thus also the branches in (4.16). On the other hand, the result 
holds for any choice of the branch. 

Proof. By a change of variables and scaling, i.e. by using the auxiliary 
function 

T-i / \ f{z + rw) - f{z) 

F{w) = — — — , weC, 
f{z + r) - f(z) 

we may assume that B = B and that / fixes the points and 1. With 

f(z) 

this scaling we are also reduced to the case where the branch of log is 
determined by the condition log/(l) = 0. 
Use then the Stoilow factorization, 



f(z) = <pof (z), 



(4.17) 



where ip is conformal on /o(2D) and /o : C — > C is a principal quasiconformal 
mapping, with fJ,f (z) = Hf{z) for \z\ < 2 and [J>f (z) = for \z\ > 2. In 
particular, 

K- 1 



(fo)z = ll{z) (f )i 



\K Z )\ ^ k X2o{z), k 



K + l 



< 1. 



Considering first the inner factor, one applies Theorem 4.3 to /o(2z)/2, 
with complex parameter 1 ^ \p — 1| ^ |. However, if we take \p— 1| = ^, at 
points where \fi(z)\ = k = ^q^j, we have p pk = 1 + pk and the integrand in 
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(4.12) vanishes, hence the estimate becomes useless for (4.16). It is for this 



reason that we need to assume the strict inequality in (|4.15). 
Setting p 



1 + 

K-l 



jr 2 ^ , we then have the strict inequalities 2 < p < 



1 + 1/k, k = f=i. Now p(z) = 1 and ( |4.12[ ) gives 



1 

7T 



9/o 



^ 4 



max{l,(l + fc) 1 " Re ^} 



(4.18) 



1 - k(p - 1) 

On the other hand, [21 (3.35)] shows for any quasisymmetric map g : Q 
O' the estimate 



\g( z ) - g(zo)\ ^ 



(20, r) 



\dg\, z G B(zo,r) C U, 



where rj(t) is the modulus of quasisymmetry as in (2.2). Further, qua- 
sisymmetry with Koebe distortion or [2, (2.61)] gives diam(/o-B) ^ c\{K). 
Therefore f B \df \^ c{K). 
Consequently, 



« K)2 < ( L 



dk 



0/2 



dk 



1-/3/2 



B 



dk 



B 



dk 



2-0 



The requirement (4.15) holds for j3 if and only if it does for 2-/3, and 



therefore (4.18) gives, too, the lower bound 

c(K,/3) ^ 



B 



dfo 



For the outer factor in the Stoilow factorization we apply Lemma 3.2 
U = / (2D). This gives 

ip(x) - <p(w) 



m 



log- 



log <p'(w) 



^ I0pa(x,w), 



x — w 

Choosing x$ = /o(l), wo = /o(0) it follows that 

tp(x ) - <p(w ) 



x,w G Q. 



log 



log/ (l)-log/o(0)| ^8(K-1) 



Xq - Wo 

where the last estimate was shown in (|3.7[). Consequently, |log (p'(wo)\ ^ 



C{K). To complete the argument note that as in the proof of Lemma 3.2 



the function | log (p'(z)\ is uniformly Lipschitz with respect to the hyperbolic 
metric of Q. Since fo(0) has hyberbolic diameter in Q bounded in terms 
of K only, | logy/(z)| ^ c(K) < 00 for z G fo(^)- With the chain rule and 
(4.17) we finally have 



B 
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Cor 4.10 



Cor 5.4 




Higher Integrability 



Figure 5. The elliptical integrability region in Theorem 4.9 
and its various consequences. 



proving the claim. 



□ 



As for higher integrability with real exponents [21 Section 13.4.1], the 
theorem can be interpreted in terms of the Muckenhoupt ^p-weights. 



Corollary 4.8. Suppose 1 < p < oo and f : 



is a K-quasiconformal 



mapping. If the exponent (3 G C satisfies both (|4.15|) and the dual condition 

(4.19) 



+ |/3 + 20-l)| <2-|^(p-l), 



K - 1 



then 



G A, 



Moreover, for each 1 < p < oo, outside this range of exponents the con- 
clusion fails for some K-quasiconformal mapping f : C — > C. 

As a last remark, the Stoilow factorization works as well for mappings 



defined in proper sub domains of C, and therefore arguing as in Theorem 4.6 



we have the following local higher integrability stated as Theorem 1.2 in the 
introduction. 

Theorem 4.9. Suppose f is a K-quasiconformal map on a domain il.cC. 
Then for any exponent (3 G C in the critical ellipse 

K + l 



+ |/3-2| < 2 



K-l 



we have 



ft 



(4.20) 
(4.21) 



Sharpness of the previous result is seen by testing with maps (3.3) 



The previous theorem includes a number of special cases which are worth 
explicating, see Figure [5} The major axis corresponds to Higher Integrabil- 
ity of pp. Other boundary points on the ellipse uncover new phenomena. 
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These will address exponential integrability of the argument for quasicon- 



formal maps (Corollary 4.10) as well for bilipschitz maps (Theorem 6.2) and 



rotational multifractal spectrum ( Corollary |5 . 4| for quasiconformal maps and 



Theorem 6.3 for bilipschitz maps). 



As a first special case with a purely imaginary exponent (5 we obtain 



(cf. Corollary 1.3) 



Corollary 4.10. Suppose f is a K -quasiconformal map on a domain ficC. 
Then 

AK 



,b\ arg /, | 



E L\ oc for all positive b < 



K 2 -l 

The result is optimal in the sense that it may fail with b = ^¥z\ f or some 



K -quasiconformal f. Such an example is provided by (3.3) with the choice 

1 



Remark 4.11. We have chosen to derive Theorem|4.9|from the precise weighted 



estimates of Theorem 4.3 However, to obtain the optimal exponent of inte- 



grability P in Theorem 4.9 one may apply the interpolation lemma in many 
different ways. For instance, we could base the argument on the standard 
holomorphic flow given by fj,\ = A • jr and the analytic family f£. In order to 
obtain uniform j£f 2 -bounds in this setting, one needs to restrict the motion 
to {|A| < 1 — e} and use quasisymmetry. 



5. Multifractal spectra 

The multifractal spectrum of a Radon measure \i on W 1 is (on the intu- 
itive level) usually defined as the Hausdorff dimension of the set of points 
x E M n for which (j,(B(x, r)) ~ r a for small radii. The rigorous definition has 
to be done carefully, and actually there are various notions of multifractal 
spectrum, see e.g. [9] or [18]. For a homeomorphism / : M n — > M. n the nat- 
ural counterpart is the multifractal spectrum of the induced push-forward 
measure [i = f*{dx). In this spirit the multifractal spectra of quasisymmet- 
ric maps of the real line was studied by the third author and Smirnov in 
|2U| . Also closely related is Binder's work [6] on the mixed integral means 
spectrum of conformal maps. 
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In this section we apply the results of the complex integrability of the 
gradient f z from the previous section, to analyze the multifractal spectra of 
a X-quasiconformal map /: C — > C. The complex integrability allows us to 
consider even the joint multifractal behaviour with respect to both rotation 
and stretching. 

According to the discussion in Section 3.2, fix a > and 7 £ M and 



consider points z £ C with the following property: 

There is a decreasing sequence (depending on z) of radii {rk)k^li 
with r/j — > 0, such that simultaneously 

\og\f(z + r k )-f(z)\ 



a 



lim 

k— ¥00 



logr fc 
arg(/(z + r k ) 



k^L log\f(z + r k 



/(*)) 

m\ 



(5.1) 



One should observe that a single point can satisfy the condition ( |5.1| ) for 
several different values of (0,7). In (5.1) the quantity \f(z + r k ) — f(z)\ 



measures stretching in the direction of the positive real axis, but by the 
quasisymmetry of / any other fixed direction gives the same result, and 
thus definition is quite robust. 

Introducing the joint rotational and stretching multifractal spectrum for 
the class of all i^-quasiconformal homeomorphisms, we are asking for a 



characterisation of the maximal size of the set of points satisfying (5.1), i.e. 
to determining the quantity 



Fk(ch, 7) := sup I dim-^(E') : (5.1) holds for every z S E, for 



some 



-fT-quasiconformal mapping / : C — >• C > , 
where dim% stands for the Hausdorff dimension. Note that, in view of 



Theorem 3.3 there are no points satisfying (5.1) unless r = a(l + 27) lies in 
the closed disk 



K + 



1 

K 



(5.3) 



The following theorem gives a complete description of the quasiconformal 
joint multifractal spectrum. 
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-r 2 








Figure 6. The joint multifractal spectrum Fx(a, 7) as a 
function of the variable a(l + 27). 

Theorem 5.1. Assume that the parameters a > 0, j £ M lie in the natural 



domain of definition of Fk, i.e. r = a(l + ^7) satisfies (5.3). In this range 
the joint multifractal spectrum equals 



F K (a,j) = (1 + q) 



~ a ) 2 ( K + !) 2 + 4i^a 2 7 2 



(5.4) 



Outside the range ( |5.3| ) we can set ^(0,7) = —00, since then the set 
corresponding to r = a(l + 17) is empty for any iC-quasiconformal map. 

Remark 5.2. It will be later useful to observe that as a function of the 



variable r = a(l + 27) the function (5.4) is determined as the unique 'cone'- 



like function on the closed disc (5.3) with the properties: the function takes 



the value 2 at r = 1, vanishes on the boundary of the disk (5.3) and is linear 



on each line segment that joins 1 to the boundary circle. Figure [6] gives an 
illustration of the graph of Fx(ct, 7). 



Proof of Theorem 5.1 We begin with upper estimates for the spectrum 
Fx(a,j). Assume that / : C — > C is X-quasiconformal and write for any 
a > and 7 G R 



Sf(a, 7) :={z£C : z satisfies (5.1) for some radii r& — > 0}. (5.5) 



Our task is to estimate dim-^(5/(a, 7)); obviously, it is enough to estimate 
the Hausdorff dimension of Sf (a, 7) n D. 

Next, we apply the complex integrability of the gradient f z established 
in the previous section. For our present purposes the most suitable form is 



given by Theorem 4.6, which in particular states for any exponent /3 in the 
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critical ellipse |/3| + \/3 — 2| < | and for any disk B(z, r) C B(0, 2) that 



f(z + r)-f(z) 



sC c(K,f3)r~ 



B 



(5.6) 



For consistency, recall that above the argument of both f z and of the dif- 
ferences f(z + r) — f(z) are obtained from a fixed choice of a branch of the 
function log [(f(w) — f(z)/{w — z)), defined for (w,z) G C 2 with w ^ z; 
the same branch is used for every z G Sf(a,'y) in the definition (5.7) below. 

Assume then that (a, 7) satisfies (5.3) with strict inequality, as it is clearly 
enough to consider this case. Fix arbitrarily small e G (0,a). By definition, 
we may select for any z G S/(a, 7) n D a radius r z G (0, e) so that 

log|/(z + r*)-/(z)| 

- G (a — e, a + e) and 

(5.7) 

G (7-e,7 + e). 



log r 2 

arg(/(z + r 2 ) - f{z)) 



log\f(z + r z )-f(z)\ 

Vitali's covering lemma allows us to select countably many points z n so that 
the discs B n := B(z n ,r n ) with r n := r Zn are disjoint and Sf(a, r y) nB C 

We fix an arbitrary f3 from the open ellipse |/3| + \j3 — 2| < § and observe 
that (5.7) together with (5.6) yields for any of the discs B(z n ,r n ) 



„2+(a-l)Re/3-a 7 Im/3+0(e) 



Above the exponent 0(e) is uniform in n. We thus obtain 



f(z n + r n ) - f{z n )y 









r 


«= / 


fS 




JB(z n ,r n ) 





y^ r 2+(a-l) Re/3-« 7 Im/3+0(e) ^ VV/3) / 
n n J £ 

<c(/3) / |(//|<(X). 



B(^ n ,r- n ) 



I (A 



(5.8) 



2D' 



Since Sf(a, r y) fll C U n 5-B n , and e > can be taken arbitrarily small, it 
immediately follows that 2 + (a — 1) Re (3 — oey Im/3 is an upper bound 
for dim-^(S , j(o:, 7)), for any /3 in the critical ellipse. As / was an arbitrary 
iT-quasiconformal map we infer that 



F K (a,j) < inf{2 + (a - 1) Re /3- cry Im/3}, 



(5.9) 
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where the supremum is taken over the set |/3| + \j3 — 2| < |. Using the 
parametrization of the ellipse in condition (4.7) we get equivalently 

( a- 1 Vl - k 2 . J 



mm 



1 + a + 



A: 



cos ( 



«7 



k 



sm ( 



(1 + a) - ^(i-a)^!!-^ 2 )"^ 2 



(1 + a) 



V(l-a) 2 (K + l) 2 + 4Ka 2 7 2 



if - 1 



ELS k 



K+l- 



In fact, as soon as we have the equality in (5.9), from this representation 
one most directly sees the cone like property of ^(0,7), as discussed in 
Remark 15.21 

It remains to find lower bounds for Fk(o(,j), and for this we provide 



examples verifying the optimality of the estimate (5.9). The examples are 



constructed by iterating the map (3.8) in a self-similar manner inside in- 
terlaced annuli that form a Cantor like structure. The quasiconformal map 
spirals only inside the annuli, elsewhere it is a similarity. The reader may 
compare the construction with that in [21 Thm. 13.6.1]). 

Again, fix K > 1 and a pair (a, 7) 6 (0, 00) x R such that r = a(l + ry) 
satisfies (5.3) with a strict inequality. This initial knowledge allows us to fix 
the auxiliary parameters t and parameters (ao,7o) via declaring t to be the 
smallest positive number such that 



t 1 (a(l + 17) - 1) + 1 =: a (l + Ho) 6 B K , 



(5.10) 



where 



B K 



T <= 



K + 



K 




Clearly t G (0, 1). There will be an additional free parameter r £ (0, 1/8) 
that we fix so small that r < 2 _4 r 4 , and denote 



„7a/ao70 



(5.11) 



We next construct a iT-quasiconformal map (f) (that depends on the pa- 
rameters ao,7o and r) by suitably iterating the model map (3.3) with pa- 
rameters (ao,7o)- For that end we select first inductively an infinite collec- 
tion of annuli, partitioned into levels j ^ 0, in such a way that the level 
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j annuli will have outer radius r J and inner radius sr 3 . In level there 
is only one annulus, i.e. 5(0,1) \ B(0,s). Assume that we have already 
constructed all the annuli up to level j — 1, where j ^ 1. For each an- 



nulus of the level j — 1, say for B(zQ,r : ' 1 ) \ B(zo,sri 1 ), we may pick 
N := (L(s/2r)J) 2 = ( [ r (^-i)/2j ) 2 ;> r 2 (*-!)/8 disjoint discs B(z m ,ri), with 
m = 1, 2, . . . , N, all lying inside the punctured disc B(zq, sr^~ l ) \ {zo}. The 
corresponding level j annuli are B(z m ,r J ) \ B(z m , sri) for m = 1,...,N. 
By performing this operation for all annuli of the level j — 1 we obtain the 
complete collection of annuli of level j, and their total number is N J . This 
collection can be written as {B \ sB : B 6 Aj}, where Aj stands for the 
set of all the outer discs of the level j annuli. Obviously this construction 
can be made via suitable similitudes so that the intersection 

E:=f]({jB) (5.12) 

j=i BeAj 

becomes a self-similar Cantor set in the plane. 

Given an arbitrary annulus B \ sB, with B = B(w,R) we define the 
corresponding rotation map ipg by setting 



» + BSS3( b ^ 1 ^ <I+ ' T0, V^ B B\,B (5.13) 
continuous similarity extension if z 6 sB. 

The quasiconformal mapping (j) is then defined via an inductive construc- 
tion, see Figure [7] for illustration. First set 4>q{z) = z and assume that 4>j-i 
is already defined. Then choose 

JK ' I ^_x(B)(0i-i(*)) lf z 6 B with B G Aj. v ' 

The new spiralling introduced via <j>j takes place in the set where 4>j-i is 
conformal (actually even complex linear). Hence, the fact ao(l + «7o) £ 
implies that each cj)j is X-quasiconformal, and so is our final map (j), where 

<j>{z) := lim (f)j(z). (5.15) 

j->oo 

Let us then consider how maps fixed level j ^ 1 balls. Any ball B £ .Aj 
has radius rj := r J and it is mapped into a ball B' = <ft(B) with radius r'- } 



where the definition of 4> together with (5.11) and (5.10) yields that 



r l, = s i(«0-l) rj = r i*K-l) = r j(a-l) = (^a ^ 5J gj 
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Figure 7. Illustration of the mapping 4> constructed in (5.15) 



Hence the stretching of <j) has the desired order with respect to the center 
point of B. Moreover, this is also true for the rotation of <j) since if B = 



B{z,r J ) we obtain directly by construction and using normalization (2.7) 
that 

arg(/(z + r°) - f(z)) = j-f a logs = -flogr' j . (5.17) 

Every point in E is inside a disc of level j. It is not quite the center point, 
but at a finite hyperbolic distance from it. Thus the above observations 
together with the robustness of the definition of pointwise rates of rotation 



and stretching (based on Theorem 3.1), allow us to conclude that 



E C S<f,(a,j). 

On the other hand, the Hausdorff dimension r := dim%(i?) is computed 
from the equation Nr T = 1, and by recalling that N = ([r^ _1 )/2j) 2 we 
obtain in the limit r — > + that r — > 2(1 — t). Hence 

F K (a,j) > 2(1 -t). 



According to definition (5.10) and the cone-type characterization observed 



in Remark 5.2, this exactly means that Fk{oi^) has the right lower bound, 



and the proof of Theorem 5.1 is complete. □ 



Combining general properties of holomorphic motions with the above qua- 
siconformal multifractal bounds quickly gives 



Proof of Theorem 1.4 Suppose f : Dx£- >■ C is a holomorphic motion of a 
set £cC. By Slodkowski's generalized A-lemma |22j ^ extends to a motion 
of the whole complex plane, while the original A-lemma of Mane, Sad and 
Sullivan [T7j proves the extended map ^x(z) = ^(A, z) to be quasiconformal 
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in C. A Schwarz- lemma type argument shows that 



1 + 



1-|A|' 



A G 



for details see e.g. [21 pp. 303-304]. The dimension bounds (1.6) hence 
follow from Theorem 15. II □ 



We next consider the upper and lower stretching exponents of / at point 
z defined by 

log|/(s + r)-/(z)| 



a f ( 



a f [z) 



lim sup 

r-»0+ 
lim inf 

r->0+ 



logr 

log\f(z + r)-f(z)\ 
log r 



In a similar manner, the upper and lower rates of rotation are given by 

arg(f(z + r) - f{z)) 



7/(z) = lim sup 



l f (z) 



log|/(* + r)-/(*)|' 

. arg(/(z + r) - f(z)) 
llmmf 1 YTi — , — \ TTV\- 

r ^o+ log \f(z + r) - f(z)\ 



(5.18) 
(5.19) 



In the following result the novelty is the estimate (5.21 ) for the Hausdorff 
dimension of the set where each point has a prescribed rotation index. This 
estimate also shows that although our method considers simultaneous ro- 
tation and stretching, it is capable of producing optimal estimates for the 
pure rotational multifractal spectra. 



Corollary 5.3. Let f : C 

Then 



K+i- 



dim^ ^{z : aj(z) = a or &f{z) = a}) 
for any a G [K~ , K]. Moreover, 



C be a K -quasiconformal map and k — lx 

1 



< 1 + a 



k 



1 



a 



dim n (iz : 7/ (z) = 7 or j f (z) = 7} J < 2 - -— = 



k 



k 



(5.20) 



(5.21) 



for any 7 with | 7 | ^ (K — i^T _1 )/2. 

If either a or 7 Zies outside the given interval, then there are no points 
with exponent a or index 7, respectively. Both estimates are optimal. 



Proof. The first statement (5.20) is deduced by a small modification of the 
proof of Theorem 



5.1 



For each z G otf(z) (resp. z G &f(z)) one chooses 
a radius r z that satisfies the first condition in (5.7) and, in the later stage 
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of the proof, employs only real parameters f3 in the allowed range (3 E 
(1 — 1/k, 1 + 1/k). In view of ( 5.8 ) and ( 5.9 ) the Hausdorff dimension of the 



left hand side set in bounded by 

inf {2 + (a-l)/3}, 
/3e(l-i/fc,i+i/fc) 

and the claim follows. That the sets in question are empty for values 
a [i-T -1 ,^] follows immediately from Theorem 3.3, and the optimality 
is obtained via considering maps (5.15) with 7 = 0. 



Towards the second statement , we first fix 7 > and write 

E := {z : 7(2:) = 7 or 7(2;) = 7}. 



(5.22) 



We claim that E C U a g/ Sj(a, 7), where / is the interval of allowed a, 
i.e. the values of a so that the pair (0,7) satisfies (5.3). Namely, if e.g. 
7(2) = 7, we may pick a sequence — > such that the second condition in 



(5.1 ) is satisfied, and by moving to a further subsequence the first condition 



(5.1) holds as well, with some allowed value for a. One observes that the 



proof of Theorem 5.1 is quite robust and yields immediately for any (a', 7) 



satisfying (5.3) and for any e > the estimate 

dim H ( [J Sf(a,i)) ^ F K (a',i) + ce 

a'—E^a^a'+e 

with a uniform constant c = 0(7, K). By covering the interval / with finitely 
many intervals (a' — e ^ a ^ a 1 + e) and since e > is arbitrary we deduce 
that 

dim. H {E) < supF^(a',7). 
a'ei 



This estimate is obviously optimal in view of the example (5.15) 



In order to compute the above supremum one may shorten computations 
by recalling from Remark 5.2 that Fx(a, 7) = 2(1 — t), with a(l + ij) — 1 = 
t(ao(l + ?7o) — 1), where ao(l + 270) is a boundary point of the disc (5.3). 



We parametrize the boundary by ao(l + ^70) = A + acos8 + iasin# where 
A := (K + l/K)/2 and a := (K - l/K)/2. Then t is determined from the 
condition 

£ao7o ta sin 6 



1 



l + i(a -l) 1 - t + t(A + a cos 1 
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and we get t = (l + (a/7) sin 9 — a cos 9 — A) 1 . The minimal value of this 
quantity is obviously 

1 



1 -A + 0^1 + 7-2 ' 
which yields the stated dimension bound. 

The emptiness of the sets under consideration in the case where the pa- 



rameters lie outside the stated ranges follows again from Theorem 3.1 Here 



one notes that the maximal slope for lines through origin that intersect the 



closed disc defined by condition (5.3) equals (K — K 1 )/2. □ 



One may also bound the size of the image of the set with prescribed 
rotation rate as follows: 

Corollary 5.4. For any K -quasiconformal map f : C — > C one has 

( \ AK 

dim w [J{z : J f {z) = 7 or j_ f {z) = 7} J < 2 - _ - | 7 | 

for any 7 with \^\ ^ (K — K~ 1 )/2. 

For 7 outside the interval there are no points with this rotation rate. The 
result is optimal. 



Proof. We again modify the proof of Theorem 5.1 Choose for each point 



z G E (where E was defined in (5.12[)) a radius r z so that the second 



condition in (5.7) holds, and pick the disjoint balls B n = B(z n ,r n ) as before. 
Define a n through r" n = \f(z n + r n ) - f(z)\ =: r' n . In ( pT8| ) and ( [5T9| ) apply 
only the values /3 = 2+iftwith |/3— 2|+|/3| < 2/k, or equivalently |6| < k^ 1 —k. 



In this situation (5.8) can be written in terms of the radii r' n in the form 
Y2n r 'rr' rb+0 ^ By quasisymmetry, the balls cB (f (z n ) , r' n ) cover the 



image f(E), where c depends only on K. Hence the analogue of (5.9) takes 
the form 

dim w (£) < inf (2 - 76) = 2 - b^AT 1 - k), 

|6|<(fc _1 -A:) 

which yields the stated estimate. 

Optimality is verified by considering the map <fi from (5.15) with the 
stretching and rotation indices (0(7), 7), with I7I < (K — K~ 1 )/2, and 
where the judiciously chosen value of a equals 

"(7) := , n ,, ■ 
1 + \^\k 
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One also notes that dim(0(£')) = d\m{E)/a as one easily computes by 
observing that 4>{E) is likewise self-similar. It is of interest to observe that, 
independently of the value of 7 > 0, the half- line determined by the points 
1 and 0(7) (1 + ij) intersects dB K at the point + jt%2- □ 

Remark 5.5. One could of course also ask for bounds of the Hausdorff di- 



mension of the image in the contexts of Theorem 5.1 or of the first part of 
Corollary |5.3[ In these cases the optimal bound is obtained by multiplying 
the already obtained results by a -1 , since one may cover the image by balls 
B(zj,cr®~ £ ), where the disks B(zj,rj) are as in the respective proofs. 

It is also reasonable to study Minkowski type multifractal spectra by 
considering dimension estimates of the type 

n / \ y ,. log N(r, a, 7, e) 
Df M (a,j) : = hmhmsup . , 

e^O r _». I log r I 

where N(r, a, 7, e) is a maximum number of disjoint disks B n = B(z n ,r) 
with center z n 6 D, such that the stretching and rotation of / on B n satisfy 



(5.7) (with r z replaced by r). The proof of Theorem 5.1 applies with mere 



cosmetic changes and yields 

Corollary 5.6. Let f : C — > C be a K -quasiconformal map. Then 



V(l-q)2(if + l)2 + 4ifaV 
Df tM {a,V ^ 1 + a - ^ — — . 



whenever r = 0(1 + ^7) lies in the disk (5.3). The estimate is the best 
possible. 

Remark 5.7. By standard localization and Stoilow factorisation, the results 
in this section generalize for /f-quasiconformal maps between arbitrary do- 
mains. 



Remark 5.8. Obviously Theorem |5.1| can be equally well understood as a 
statement on the combined ordinary multifractal spectrum of the pull-back 
measure fi := f*(dx) and the rotation spectrum. 
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6. BlLIPSCHITZ MAPS AND ROTATION 

In this section, we apply our quasiconformal methods to bilipschitz maps. 
These have, by the very definition, trivial stretching but may exhibit non- 
trivial rotation. On the other hand, extremal quasiconformal mappings for 
various "stretching problems" often exhibit no rotation at all. This leads 
to a number of dual analogies, as described in Table [I] in the introduction. 
Below we uncover the various entries from this dictionary. 

As discussed in Section [2. 1[ orientation preserving L-bilipschitz mappings 
/: f2 — > C are L 2 -quasiconformal (and if / is not orientation preserving, 
then its conjugate / is). An archetypical example of a bilipschitz map with 
spiralling behaviour is the logarithmic spiral map 

s y (z) = z\zP, 7 £l, (6.1) 

that we already discussed before. Indeed, it is not hard to verify that 
s 7 : D — > D is L-bilipschitz where L ^ 1 satisfies 

i-i-H- 

In the following, we study rotational behaviour of general planar bilipschitz 
mappings in various points of views. All of our results will be consequences 
of the quasiconformal theory. Nevertheless, in this way we obtain sharp 
results even in the bilipschitz category. One way to explain this phenomenon 
is that often the maps extremal in the bilipschitz category, such as the 
logarithmic spiral map (and its iterated variants), are not only bilipschitz 
but also area-preserving; thus there is an exact correspondence between the 
optimal quasiconformal and bilipschitz constants of the form K = I? . 

6.1. John's problem. Denote by A = A(r, R) = {z £ C : r < \z\ < R} the 
annulus with radii < r < R < oo. Let /: A — > A be an L-bilipschitz map 
identity f(z) = z on the outer boundary \z\ = R, with a prescribed rotation 
(parametrized by 7 £ 1) on the inner boundary, 

f(z) = «*rk>g(fl/r) ) \ z \=r. 

John's problem asks for constraints on the bilipschitz constant L for such 
a map to exist. Note that the analogous problem of finding quasiconfor- 
mal maps between annuli (of different conformal modulus) is attributed to 
Grotzsch. In his work |15j F. John obtained some quantitative results in 
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the asymptotic regime L — > 1, while the complete solution was given by 
Gutlyanskh and Martio [12] : the rotation parameter 7 and the bilipschitz 
constant L need to satisfy 

| 7 |^L-i (6.2) 



As discussed above, the spiral maps such as in (6.1 ) provide extremal exam- 
ples. For related recent results in the class of mappings of finite distortion, 
see [1]. 



We begin by generalizing (6.2) to all L-bilipschitz maps between arbitrary 
domains. 

6.2. Pointwise rotation. Recall the upper and lower rates of rotation from 
Section [5| definitions (5.18). For bilipschitz mappings it is equivalent to use 



the following convenient formulations 

arg(/(z + r) - f(z)) 
logr 

arg(/(z + r) - f(z)) 



'jf(z) = limsup 
1 r^0+ l °g r 



7 (z) = liminf 
-J r^o+ log r 



Proposition 6.1. Let f'.Q—tQ! be an L-bilipschitz homeomorphism be- 
tween planar domains. Then the rates of rotation satisfy the following point- 
wise bound 

1 

T 



17/(3)1,17/2)1 < L 



Z £ ft. 



The spiral map (6.1) shows that this is best possible in general. 



Proof. After taking the conjugate if necessary, the L-bilipschitz map be- 
comes L 2 -quasiconformal. We observe that the bilipschitz property of / 



ensures that the limit (3.11) is equal to one at every point regardless of the 



subsequence of radii chosen. Hence an application of Theorem 3.3 yields for 
(any) rate of rotation 7 the inequality 

1 



1-1 



K 



1 
K 



< 



K 



1 
K 



As K ^ L 2 , a simplification yields the stated bound. 



□ 
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The slightly weaker estimate I7I ^ VL 2 — 1 was derived by Freedman and 
He in connection with the bilipschitz factoring problem for the map s~. 



We will see in Section 6.5 how the optimal form of Theorem 6.1 leads to an 



optimal answer of the factoring problem, as well. 



Theorem 6.1 is sharp as a pointwise estimate, but one expects that ex- 
tremal spiralling behaviour cannot simultaneously occur at many places. 
This is the theme we discuss next. 

6.3. BMO and Exponential integrability. On a historical side, the func- 
tion space of bounded mean oscillation was originally introduced by John 
and Nirenberg |15[ [16] exactly in the context of rotational phenomena of 
bilipschitz maps, i.e. John's problem above. Roughly speaking, John estab- 



lished a discrete variant of the fact, c.f. Proposition 3.8, that 

arg/ z e BMO. 

This coupled with the John-Nirenberg lemma leads then to bounds on the 



rotation problem considered in section 6.1 

In view of the previous discussion it is natural to ask for the best possi- 
ble exponential integrability of the argument arg/ 2 . As an application of 



Theorem 4.9 we obtain immediately 



Theorem 6.2. Let f: f2 — )■ C be an L-bilipschitz mapping. Then for any 
^ b < jjrzi we have 

exp(&|arg/,|)e£L(fi). (6.3) 

Furthermore, the integrability fails at the borderline b = f or some L- 
bilipschitz mapping f . 

Proof. Since \f z \ is uniformly bounded above and below, only the imaginary 



part of f3 plays now a role in Theorem 4.9 But for /3's inside the critical 



ellipse |/3 1 + \/3 — 2| < 2 ^ 2 ^j , the supremum of the imaginary part equals 



A_l (and infimum equals — J L t ). The optimality is seen by considering 



the model map (6.1) with L — 1/L = □ 



6.4. Multifractal spectrum. Interpreting Theorem 6.2 more geometri- 



cally we use it for multifractal bounds on rotation. Indeed, by specializing 



Theorem 5.1 to the bilipschitz case, that is, by setting a = 1 and K = L 2 , 



we obtain the following characterisation 
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T 2 



1/K 



- 

K -(L-1/L) 



L-1/L 



Figure 8. Multifractal spectra: quasiconformal stretching 
(left) and bilipschitz rotation (right). 



Theorem 6.3. Let f : 0, — > f2' be an L-bilipschitz mapping, where and f2' 
are planar domains. Then we have the optimal bounds 

2L 

dim H ({z£Q: j f (z)=^y or 7 (z) = j}) ^ 2 



L 2 -l 



l7l 



for every admissible I7I ^ L 



1 

L ■ 



Theorem |1.1| is an immediate corollary. The optimality of the bounds, 
i.e. that they hold as an equality for some L-bilipschitz mapping, follows 
from the construction in Theorem 5.1 with a = 1. Figure [8] above contrasts 



Corollary 5.3 and (5.20) with Theorem 6.3 Observe also that for bilipschitz 



maps the dimension of any set is preserved, whence the above estimate is 
valid also for the image of the set of prescribed rate of rotation. 

6.5. Factoring the logarithmic spiral. A basic open question in the 
study of bilipschitz mappings (in R n ) is whether such a map can be rep- 
resented as a composition of (1 + e)-bilipschitz mappings, for any e > 0. 
The factoring is known only in dimension n = 1; see |10j for recent general 
results on this theme. Towards this problem Freedman and He studied 



the factoring of the logarithmic spiral map (6.1). As an application of our 



pointwise rotation estimates we will revisit their question. 

One way to factor s 7 to maps of small bilipschitz distortion is to simply 
write it as a composition of slower spirals: 



S'y — ^""yo 



*70 



o s 



'70 ) 



N terms 



where 70 = 7/-ZV. The next Theorem says that this is the most efficient way. 
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Theorem 6.4. Let s 7 : D — > D 6e factored as s 7 = /at o /^v-i ° • • • /i, where 
each fi is an Lo-bilipschitz map of a closed Jordan domain in M? , Lq > 1. 



Then the number of factors needed is at least N ^ 



h\ 

Lq — j 



Freedman and He proved the same result with lower bound \ j\/ \/L^ — 1, 
and the improvement above was observed in [12] for factoring within a special 
class of bilipschitz maps. 

Proof. We follow [IT], where the only adjustment needed is the sharp form 



of Theorem 6.1 For simplicity, we assume 7 > 0, the 7 < case being 
similar. The crucial observation of [11] is the subadditivity of the rate of 
rotation under composition. In our notation, for / and g bilipschitz maps, 

J 9 of(z) < 7f(z) + j g (f(z)), 

provided that 7 y go j-(z) > 0. Since 7 S7 (0) = 7, the repeated application of 
subadditivity and the estimate of Theorem | 6.1| implies 

1 



7 ^ N L G 

as required. □ 



Remark 6.5. The content of Theorem 6.4 is that bilipschitz factoring even 
when exists might need exponentially more factors than optimal quasicon- 
formal factoring. One may visualize the difference in this particular example 
by considering the maps f T {z) = A where r = a(l + ij), with the 
parameter space H = {Rer > 0}. Here the logarithm of the quasiconfor- 
mal distortion logiT(/ r o Z^ 1 ) equals the hyperbolic distance d^{r,r'). If 
we now start at s 7 (z) = z\z\ 11 , travelling to the identity along the hyper- 
bolic geodesic in H provides us the optimal quasiconformal factoring, while 
bilipschitz factoring requires travelling along the horocycle a = 1. 
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